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Abstract 

A method is presented for proving upper bounds on the moments of the position op- 
erator when the dynamics of quantum wavepackets is governed by a random (possibly 
C^l ! correlated) Jacobi matrix. As an application, one obtains sharp upper bounds on the dif- 

fusion exponents for random polymer models, coinciding with the lower bounds obtained 
in a prior work. The second application is an elementary argument (not using multiscale 
analysis or the Aizenman-Molchanov method) showing that under the condition of uni- 
formly positive Lyapunov exponents, the moments of the position operator grow at most 
logarithmically in time. 



o 

1 Introduction 



One of the fundamental questions of quantum mechanics concerns the spreading of an initially 
localized wave packet under the time evolution e~ ltH associated to a Schrodinger operator 
H. If the physical space is M. d or Z d and the position operator is denoted by X, the spreading 
can be quantified using the time-averaged moments of X (or equivalently the moments of the 
associated classical probability distribution): 

/OO J. 
^e-%|e^|Xre-^|0>, g>0. (1) 

It is well known that for short-range operators H, the moments cannot grow faster than ballis- 
tically, that is Mj, < C(q)T q . The growth actually is ballistic in typical scattering situations 
and for periodic operators H describing Bloch electrons. On the other hand, if the moments 
M|. are bounded uniformly in time, one speaks of dynamical localization. This can be proven in 
the regime of Anderson localization for random operators, but also for certain almost-periodic 
operators (see [Jitj for a review). There are many models where the moments Mj, exhibit 
some non-trivial power law behavior. If Mj, ~ T q ^ 2 , the quantum motion is diffusive, and any 
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other asymptotic growth behavior is called anomalous diffusion. In order to distinguish various 
anomalous diffusive motions, one defines the diffusion exponents 

8+ = li msup lQ g( M ^ 3- = hminf log(M ^ a > (2) 

P * T?~ P log(T«) ' P * ™ log(T9) ' q [ ] 

If the limit exists, we write (3 q = (3^ = (3~ . The diffusion exponents correspond to the Levy- 
Khinchin classification of Levy flights in classical probability, however, we stress that the quan- 
tum anomalous diffusion does not result from a probabilistic dynamics, but rather from a 
Hamiltonian one. It is due to delicate quantum interference phenomena. The ballistic bound 
implies < (3^ < 1 and convexity inequalities show that f3f- is non-decreasing in q. In the 
regime of dynamical localization (3 q = and for quantum diffusion f3 q — |. Anomalous dif- 
fusion corresponds to all other values of (3^. Typically (3^ is then also varying with q which 
reflects a rich multiscale behavior of the wave packet spreading. Such anomalous diffusion was 
exhibited numerically in several almost-periodic Jacobi matrices having singular continuous 
spectra (e.g. Fibonacci and critical Harper operator), and also some random and sparse Jacobi 
matrices. 

It is a challenging problem of mathematical physics to calculate the diffusion exponents for 
a given Schrodinger operator, in particular, when the quantum motion is anomalous diffusive. 
In this work we accomplish this for the so-called random polymer models, a random Jacobi 
matrix described in detail in the next section, and show that 

(3 q = max jo, 1-^} ?>0, (3) 

see Theorem [31 This result confirms the heuristics and numerical results of Dunlap, Wu and 
Phillips [D WP] for the random dimer model, the prototype of a the random polymer. The latter 
model was introduced and analyzed in our prior work in collaboration with G. Stolz [JSSJ, which 
already contained a rigorous proof of the lower bound (3~ > 1 — i. In this work we hence focus 
on the upper bound, which amounts to proving quantitative localization estimates. 

Next let us discuss this result in the context of prior rigorous work on other one- dimensional 
models (Jacobi matrices) exhibiting anomalous diffusion. First of all, the Guarneri bound 
|Gua] and its subsequent improvements |Coml ILas| IGS2| IGS3| IKL| IBGT] allow to estimate the 
diffusion exponents from below in terms of various fractal dimensions of the spectral measure. 
However, those results do not allow to prove the lower bound in ((3D because, as was shown 
by de Bievre and Germinet [BGj . the random dimer model has pure-point spectrum so that 
the Hausdorff dimension vanishes and the Guarneri bound is empty. In fact, the argument in 
[ JSS] is based on a large deviation estimate on the localization length of the eigenstates near 
the so-called critical energies at which the Lyapunov exponent vanishes. 

Upper bounds on anomalous quantum diffusion were first proven for Jacobi matrices with 
self-similar spectra \GS1\ [BS] . and these bounds are even optimal for the so-called Julia matri- 
ces. Kiselev, Killip and Last |KKLj presented a technique based on subordinacy theory allowing 
to control the spread of a certain portion of the wave packet (not the fastest one and there- 
fore not the moments), and applied it to the Fibonacci model. Tcheremchantsev |Tch] proved 



2 



tight upper bounds on growing sparse potential Hamiltonians introduced in |JLj and further 
analyzed in [CM] . Recently, Damanik and Tcheremchantsev [DT] developed a transfer matrix 
based method that allows to prove upper bounds on the diffusion exponents and also applied 
it to the Fibonacci model. Another way to achieve upper bounds on the diffusion exponents in 
terms of properties of the finite size approximants (Thouless widths and eigenvalue clustering) 
was recently developed and applied to the Fibonacci operator by Breuer, Last, and Strauss 
[BLS]. In the models considered in |GS1[ IBS[ IDT[ IBLS[ ITch] the anomalous diffusion is closely 
linked to dimensional properties of the spectral measures, even though the (generalized) eigen- 
functions have to be controlled as well. As a result, the transport slows down as the fractal 
dimension of the spectral measure decreases. The origin of the anomalous transport in the 
random polymer model is of different nature. In fact, in the random polymer model only a few, 
but very extended localized states near the critical energies lead to the growth of the moments. 
Hence this model illustrates that spectral theory may be of little use for the calculation of the 
diffusion exponents. This statement is even more true if the dimension of physical space is 
higher. There are examples of three-dimensional operators with absolutely continuous spectral 
measures, but subdiffusive quantum diffusion with diffusion exponents as low as | |BeS] . 

The strategy for proving upper bounds advocated in [DTJ appears to be more efficient in the 
present context than prior techniques [GS1, KKLJ. We refine and generalize the relevant part in 
Section [3], see in particular Proposition [2j It allows to give a rather simple proof of our second 
result, namely Theorem [H which establishes a logarithmic bound Mj, < log(T) /3g , (3 > 2, under 
the condition of uniformly positive Lyapunov exponent. This result is neither new nor fully 
optimal. However, in the generality we have (any condition on the distribution of randomness, 
including e.g. Bernoulli) the Aizenman-Molchanov method [AMj cannot be applied, and the 
only technique previously available to obtain this statement was the multi-scale analysis of 
[CKM] (see also |BG[ IDSS] ). Thus our proof is significantly simpler. Moreover, one can 
argue that it captures the physically relevant effect of localization. Indeed, it was shown by 
Gordon |Gorj and del Rio, Makarov and Simon [DMSJ that a generic rank one perturbation of a 
model in the regime of strict dynamical localization (M|. bounded) leads to singular continuous 
spectrum, and therefore, by the RAGE theorem, growth of the moments. However, it was 
shown in |DJLSj that this growth can be at most logarithmic, just as proven in Theorem [TJ 
The proof of Theorem [1] constitutes essentially a part of the proof of Theorem [3j 

In the next section we present our models and results with technical details. Section [3] 
contains a general (non-random) strategy for proving the upper bounds. Section H] provides 
the proof of Theorem [1] as well as some statements that are used in Section [61 In Section [5] we 
obtain probabilistic bounds on the transfer matrices near a critical energy based on the large 
deviation estimate of [JSSJ. Section [6] contains the proof of Theorem [31 that is the identity 

Acknowledgment: This work would have been impossible without |JSSj . We are very thank- 
ful to G. Stolz for this collaboration. We also thank the anonymous referees for comments that 
improved the paper. The work of S. J. was supported in part by the NSF, grant DMS-0300974, 
and Grant No. 2002068 from the United States-Israel Binational Science Foundation (BSF), 
Jerusalem, Israel. H. S.-B. acknowledges support by the DFG. 
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2 Models and results 



A Jacobi matrix is an operator H u on £ 2 (Z) associated to the data uj = (t n ,v n ) nE %, of positive 
numbers t n and real numbers v n which we suppose to be both bounded by a constant C, and 
t n bounded away from 0. Using the Dirac notation |n) for the canonical basis in £ 2 (Z), is 
given by 

Huj \n) = t n+ i \n + 1) + v n \n) + t n \n - 1) . (4) 

Each uj is called a configuration. The set of all configurations is contained in Q = ([— C, C] x2 ) xZ . 
The left shift S is naturally defined on Q. A stochastic Jacobi matrix is a family (H^uj^q, of 
Jacobi matrices drawn according to a probability measure P on O which is invariant and 
ergodic w.r.t. to S. Furthermore, we speak of a random Jacobi matrix if P has at most finite 
distance correlations, namely there exists a finite correlation length LgN such that (t n , v n ) and 
(tm,v m ) are independent whenever \n — m\ > L. Most prominent example of a random Jacobi 
matrix is the one-dimensional Anderson model for which t n = 1 and the v n are independent 
and identically distributed so that L = 1. Random polymer models as studied in [DWPl FJSSj 
and described in more detail below provide an example of a random Jacobi matrix with finite 
distance correlations (the second crucial feature of these models is that the (t n ,v n ) only take a 
finite number of values). 

We will consider the disorder and time averaged moments of the position operator X on 
£ 2 (Z), denoted by Mf, as in ©: 

M« = ^-e-^E (0| e lH ^ \X\ q e'^ |0) , q > . (5) 
Jo T/2 

Here E denotes the average over u w.r.t. P. (Note that upper bounds on the expectation w.r.t. 
P yield upper bounds almost surely.) One may replace |0) by any other localized initial state 
(at least in ^(Z)). 

As discussed in the introduction, it is well-known that the Anderson model exhibits dynam- 
ical localization, that is M|. < C(q) < oo uniformly in T for all q > 0. A by-product of our 
analysis of the random polymer model discussed below is a simple proof of the weaker result 
that Mj, grows at most logarithmically in T whenever the Lyapunov exponent is strictly posi- 
tive. In order to define the latter, let us introduce as usual in the analysis of one- dimensional 
systems the transfer matrices at a complex energy z by 

V(n,m) =TZ_ 1 -...-T', n>m, T* = ( "*» 

Furthermore Tj(n,m) = T*{m,n)~ x for n < m and Tj(n,n) = 1. Then the Lyapunov 
exponent is 

7 (*) = lim log (K(N, 0)||) . 

N— >oo iv 

Because the t n and v n are uniformly bounded, one shows by estimating the norm of a product 
of matrices by the product of their norms that for every bounded set U C C 

\\r:(n,m)\\ < e^" m ' zeU, (6) 
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where the 71 depends on U. This implies that 7(2) < 71 for z G U . For a random Jacobi 
matrix, uniform lower bounds on 7(2) can be proven by the Furstenberg theorem (e.g. |PF] ) . 
This applies in particular to the Anderson model (also with a Bernoulli potential), and more 
generally to random Jacobi matrices with correlation length equal to 1. 

Theorem 1 Consider a random Jacobi matrix. Let the (non-random) spectrum be o~(H) C 
(£"0, Ei) Suppose that the Lyapunov exponent is strictly positive: 

l(z) > 7o > 0, ze(E ,Ei) . (7) 

Then for any (3 > 2 there exists a constant C((3, q) such that 

M| < (logT)^ + C(f3,q) . (8) 

Another self-averaging quantity associated to to ergodic Jacobi matrices is the integrated 
density of states (IDS) M(E) which can be defined by 

J M(dE) f(E) = E (0|/(# w )|0) , / g C (R) • 

The Thouless formula (e.g. |PF] ) links the Lyapunov exponent to the IDS. It implies that, if 
(J7|) holds for real E G (E , Ex), then one also has 7(2) > 70 for all z with ^te(z) G (E , Ex). 

As shown by Theorem [1] and the remark just before it, it is necessary for a random Jacobi 
matrix to have a correlation length larger than 1 for the moments My. to grow faster than 
logarithmically, so that the diffusion exponents f3 q do not vanish. That this actually happens 
for the random dimer model was discovered by Dunlap, Wu and Phillips |DWP] . 

Next let us describe in more detail the more general random polymer model considered in 
[JSS]. Given are two finite sequences t± = (t±(0), . . . , t±(L± — l)) and v± = (v±(0), . . . , v±(L± — 
1)) of real numbers, satisfying t±(l) > for all I — 0, . . . , L± — 1, L± > 1. The associated 
random polymer model is the random Jacobi matrix constructed by random juxtaposition of 
these sequences and randomizing the origin. More precisely, configurations to G Q can be 
identified with the data of a sequence of signs (a n ) n( zi and an integer < / > L ai — 1, via 
the correspondence (t n ) = (. . .,i ai (l), . . . , t ai (L ai - l),t a2 (0), . . . ,i a2 (L a2 - l),t ff3 (0), . . .) and 
similarly for (v n ), with choice of origin t = t ai (l) and vq = v ai (l). The shift is as usual, and 
the probability P is the Bernoulli measure with probabilities p + and p_ = 1 — p + combined 
with a randomization for / (cf. [JSSJ for details). The correlation length in this model is 
L = max{L + , It is now natural and convenient to consider the polymer transfer matrices 

Ti = T*. ut (T ....... T. z rm f (n , , where T* = ( ^ ~ ^ * ' ~* j ■ (9) 

± v±{L±-l),t±(L±-l) v±(0),t±(0) ' v,t I f-i- Q I V / 

Definition 1 An energy E c 6l is called critical for the random polymer model (H u ) u ^q if the 
polymer transfer matrices T± c are elliptic (i.e. \Tr(T± c )\ < 2) or equal to ±1 and commute 

[Tf c ,T_f c ] = . (10) 
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If L± = 1, the model reduces to the Bernoulli- Anderson model and there are no critical 
energies. The most studied [DWPt IBovl IBG] example is the random dimer model for which 
L + = L_ = 2 and t±(0) = t±(l) = 1, v+(0) = v+(l) = A and u_(0) = = -A for some 

A G R. This model has two critical energies E c = A and E c = — A as long as A < 1. For 
further examples we refer to [JSSj . It follows from the definition that a simultaneous change of 
coordinates reduces both T + and T_ to rotations by angles that are denoted by r\ + and r\_. 

It is immediate from (jlOjl that the Lyapunov exponent vanishes at a critical energy. Because 
the transfer matrices T± are analytic in z, it follows that there is a constant Cq such that for all 
e G C with |e| < e one has for n, m G N, 

||7^ +e (n,m)|| < e Co|e||n - m| . (11) 

In particular, \^{E c + e)\ < cq |e|. However, the correct asymptotics for the Lyapunov exponent 
is 7(£ , c + e) = 0(e 2 ). This was first shown (non-rigorously) by Bovier for the case of the random 
dimer model [Bov] . but heuristics were already given in jDWPj . The rigorous result about the 
Lyapunov exponent and also the integrated density of states M are combined in the following 
theorem. 

Theorem 2 [JSSj Suppose that E(e 2M?CT ) ^ 1 and E(e 4M?<T ) ^ 1. Then for e G R and some 
D > 0, the Lyapunov exponent of a random polymer model satisfies 

7 (£ c + e) = De 2 + 0(e 3 ) , (12) 

in the vicinity of a critical energy E c . //|| [T5 +e ,Tf c+e ] || > Ce for some C > and small e, 
one has D > 0. Moreover, the IDS M satisfies 

Af(E c + e) - Af(E c - e) = D'e + C(e 2 ) , (13) 

/or some constant D' > 0. 

Furthermore [JSSj contains explicit formulas for D and £)'. The bound D > is not 
explicitely contained in [JSSJ, but can be efficiently checked using Proposition 1 in [SSSJ. 
The statement about the IDS only requires E(e 2ir?CT ) ^ 1. Let us remark that the hypothesis 
E(e 4 * r?CT ) 7^ 1 does not hold, for example, in the special case of a random dimer model if 
A = 1/ v2- In this situation, one is confronted with an anomaly. Nevertheless, the asymptotics 
is as in f)12p and again one can calculate D explicitely [SchJ. As we did not perform the large 
deviation analysis of [JSSJ in the case of an anomaly, we retain the hypothesis of Theorem [2] 
throughout. 

The following is the main result of this work. 

Theorem 3 Suppose that E(e 2 * T?CT ) ^ 1 and E(e 4i, ' CT ) ^ 1, and that the random polymer model 
has a critical energy at which (|T2|) holds with D > 0. Then for q > 0, 

(3 q = max I , 1 - — - 
I 2g 

As already pointed out, the lower bound (3~ > 1 — j- was already proven in [JSSj . 




6 



3 A strategy for proving upper bounds on dynamics 

In this section it is not necessary for the Jacobi matrix to be random or ergodic; hence the 
index u> and the average E are suppressed. Let the notation for the Green's function be 

G z (n,m) = (n\ — - — \m) , (14) 
H — z 

where n, m G Z and z G C is not in the spectrum a(H) of H. The starting point of the analysis 
is to express the time averaged moments flSJ) in terms of the Green's function 

K = EM* / ^\G E+ H0,n)f. (15) 



n€Z 



Using the spectral theorem, this well known identity can be checked immediately by a contour 
integration. In order to decompose the expression on the r.h.s., let us introduce for < a < «i 
and E < Ei, 

. r E i f^vp 

Mr°' ai (^o,^i) = £ H' ^\G E ^(0,n)\ 2 , (16) 

T a O<\n\<T a l E ° 

and Mjl ' a (Eo,Ei) is defined similarly with the sum running over < \n\ < T a , a > 0. The 
following result also holds for higher dimensional models. 

Proposition 1 Suppose o~(H) C (E ,Ei) and a > 1. Set A = dist({E ,Ei},a(H)). Then 
there exists a constant C\ = Ci(a, A,q) such that 

|M| - Mf> a {E Q ,E l ) \ < ^. 

The proof is based on the following Combes-Thomas estimate. Even though standard, its 
proof is sufficiently short and beautiful to reproduce it. 

Lemma 1 Let A(z) = dist(z, c(H)) and C2 = (4||t|| 00 ) _1 where ||t||oo = su P n gz^- Then 

2 

\G z (n, m) \ < exp ( — arcsinh(C2 A(^)) |n — m|) . 

l\[z) 

Proof. For 77 G M, set H v = e vX He~ vX . A short calculation shows that 

\\H„ -H II < lltlL le" - e~ v \ . 



Hn-z)-^ < ( \\(H- z^W 1 -\\H V -H 



Hence one has 



Since \\(H — z) 1 \\ < A(z) 1 , the choice \r]\ = arcsinh (A(z)/(4 ||t||oo)) hence gives \\(H V — 
z)- 1 ]] <2/A(z). The bound now follows from (n\(H - z)- l \m) = e ^ m ~ n ) {n\(H n - z)- l \m). □ 

The following estimate will be used not only for the proof of Proposition [U but at several 
reprises below. 
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Lemma 2 Let A, a > 0, q > and N £ N. Let p = [^— ] where [ ■ } stands for the integer part. 
Then 



n>N 



Proof. Bounding the sum by the integral on each interval of monotonicity and then extending 
this integral to the entire range [N, oo), we obtain 



poo poo -AN a P | 

Y n*e- Ana < 2 / dxxie~ Axa < - / dy y p e~ Ay = V V - (NAV . 

Bounding the sum over j by p\ (NA + l) p completes the proof. □ 

Proof of Proposition [TJ We first consider the energies above the spectrum and set A = 
dist(-Bi, a(H)). Due to the previous two lemmata and if p = [q+ 1], 



~T S H 9 T-g exp(-arcsinh(C 2 (A + £))M) 
71 M>i 

^ / ^A^T/F arcsinh^A + S))-^ 1 ) . 
Jo 7rT A + E 

Since arcsinh(?/) > ln(y) for large y, this shows that M|,' ' oo (i? 1 , oo) < C /T for some constant 
C = C(g, A). A similar bound holds for M^ 0,oo (~ oo, _E ). Now using the imaginary part of the 
energy in Lemma [Hand the bound arcsinh(?/) > y for sufficiently small y > 0, we obtain 

M^ a '°°(E , Ei) < - [ 1 dE y \n\ q exp(-arcsinh(C 2 /T)|n|) 

71 Je |n|>T" 

For a > 1, this decreases faster than any power of T. Combining these estimates implies the 
proposition. □ 

According to Proposition [TJ and equation ( fT6l) . one now needs a good bound on the decay 
(in n) of the Green's function for complex energies in the vicinity of the spectrum. As shown 
by Damanik and Tcheremchantsev [DTj . such bounds can be obtained for Jacobi matrices in 
a very efficient way in terms of the transfer matrices. Here we give a streamlined proof of this 
statement which also works for arbitrary Jacobi matrices (the kinetic part is not necessarily 
the discrete Lapalacian) and does not contain energy dependent constants as in [DTj . It will 
be convenient to first consider such bounds for the half-line problem, which is operator (jlj) on 
£ 2 (N) with Dirichlet boundary conditions. This operator and its Green's function are denoted 
by H and G z (n,m). 
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Proposition 2 Set r = maxlJtH 2 ^ , 1 , \\t 1 \\l Q } and z = E + ^. One has the bounds 

E4 r 3 T 4 

l^(0,n)| 2 < r^r> 

and, for T > 1, 

^ N ~ max < |n| < w || T*(n, 0)||' 

Proof. Let T N be the decoupling operator at N denned by T N = t N+1 (\N)(N+l\ + \N+l)(N\) 
and set G Z N = (H — z — T N )~ l and G z = (H — z) -1 . The resolvent identity reads 

P^Z fiZ fiz "p /^z 

Thus, with the notation G z N (n,m) = (n\G z N \m), 
J2\G z (0,n)\ 2 = J2\ G N T NG z (0,n)\ 2 

n>N n>N 

= \G z N (0,N)t N+1 \ 2 J2\G z (N+l,n)\ 2 < ||t||^T 2 |G^(0, N)\ 2 , 

n>N 

since T^ 1 = ^sm(z) and \\G Z \\ < T. As the l.h.s. is decreasing in N, we therefore have 

J2\G z (0,n)\ 2 < \\t\llT 2 min K(0,n)| 2 . (17) 



n>N 



Now let Iljv = ^2n=o \ n )( n \ ^ e the projection on the states on the first N + 1 sites and set 
Hjsr = HnHII n . Then G z N (n,m), < n, m < N, are the matrix elements of the inverse of an 
(N + 1) x (N + 1) matrix H^ — z. The matrix elements are closely linked to the transfer matrix 



since a, 6, c, <i when multiplied by the n^=i ^, are the determinants of certain minors of z — Hm- 
Namely by Cramer's rule (or, alternatively, by the Stieltjes continued fraction expansion and 
geometric resolvent identity) the following identities hold for z a{H^): 

^(0,0) = i - , G Z N (N,N) = -- , 0^(0,0) = i - , 

and 

G z N (0,N) = - ii, (7^(0, iV-1) 1 ; 



to o to tjv c 

Therefore |6| < to T l a l> l c l < T M and Ml < ^o T l c l- As the matrix norm is bounded by the 
Hilbert-Schmidt norm, it follows that 

AT 2 r 2 



|T 2 (A^ + 1,0)|| 2 < 



min{|G^(0, AO| 2 ,|GVi(0, iV-l)| 2 } 
9 



By (I17p this proves the first inequality. The second one follows from the first one (coupled with 
the same statement for the left half-line) by observing that the resolvent identity gives 

G z {0,n) = G z {0,n) - G z {0, -1) t G z {0, n) . 

Therefore |G 2 (0,n)| < (1 + Tto) \G z (0,n)\ which implies the second bound. □ 



4 Logarithmic bounds in the localization phase 

In this section we provide the proof of Theorem [1] and hence suppose throughout that the 
stated hypothesis hold. The main idea is to use the given positivity of the Lyapunov exponent 
d?]), combine it with the given uniform upper bound §6§ in order to deduce good probabilistic 
estimates on the growth of the transfer matrices. This growth in turn allows to bound the 
Green's function due to Proposition [2] which then readily leads to the logarithmic upper bound 
on the moments. 

Let us set U = {z eC\E < Ue{z) < E 1 , \Qm(z)\ < 1 }. 
Lemma 3 For z eU and iVeN, the set 

n N (z) = {wGfl I ||TJ(iV,0)|| 2 > e^ N } 

satisfies 

P(<>v(:)) > 



2 7i - 7q 



Proof. Let us set P = P(Qn(z)). Due to (CO), the subadditivity of the transfer- matrix cocycle 
and the bound §6§, it follows that 

70 < log(||TJ(iV,0)||) < (l-P)^ 7 o + P-Y1, 

with 71 defined by ([6]) using U as above. This directly implies the result. □ 
Lemma 4 Let z G U and N G N. Then there is a constant C3 = 6*3(70,71) such that the set 



VLn{z) = <{ to E £1 

satisfies 



max ||T 2 (n,0)|| 2 > e C; 

0<n<N " ^ V ' y " - 



P(Cl N (z)) > 1 - e~° iN ^ . 

Proof. Let us split iV into pieces of length iVo (here and in the sections below, we suppose 
without giving further details that there is an integer number of pieces and that the boundary 
terms are treated separately). By the stationarity, on each piece [jN + 1, (j + l)iV" ), Lemma [3] 
with N = N applies. As the pieces are independent, we deduce 



uj e Q 



max ||TJ((j + l)iVo,jiVo + l)|| 2 < e^ < (l-p )^ 

0<j<N/N o 
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wherepo = 7o/(2 7i-7o)- Furthermore Tj((j + l)N ,jN + l) = Tj((j + 1)N , 0)% z (jN o , 0)" 1 . 
As A = BC implies either ||£?|| > ||^4|| 5 or ||C|| > \\A\\^ for arbitrary matrices, and = 
1 1 A 1 1 for A G SL(2, C), it therefore follows that 



P N co G ft 



max ||TJ(jiVo,0)|| 2 > e^° N ° \) > 1- (l-po)$ 
a<j<N/N ) J 

Choosing AT = ciV^ with adequate c concludes the proof. □ 

Since the above lemma applies equally well to negative integers N, it is a direct corollary 
of Lemma H] and Proposition [2] that, for T sufficiently large, 

E |G:(0,n)| 2 < 32r 4 T 6 e- CaNh , z = £ + 1 G C/ . (18) 



|n|>7V 



The following lemma, holding for arbitrary ergodic families of Jacobi matrices, is useful for 
bounding a sum similar to the one in ( fTHj) . but over |n| < N. Let B^(z) = J ji{dE) (z — E)~ l 
denote the Borel transform of a measure /x. 



Lemma 5 For any F G K, T > and N > 1, one has 

J2 ^E|G:(0,n)| 2 < — QmB„(z), z = E+-. (19) 

z — ' 71 1 71 1 



0<\n\<N 

Furthermore, for any Eq < E\ , 

/•Pi 7 I", 



£ K / ^E|G^(0,n)| 2 < TV*. (20) 

0<\n\<N 



Proof. One has 



V ^-E\G z JO,n)\ 2 < N qy y -i=E(0|- ^ r |n) (n|- ?- -|0) 

0<|n|<V n&L T w T 

< iV 9 E — (0|- — r |0) , 

which by the spectral theorem gives (|T9|) . Using 

El -3m^(£ + l r 1 ) < [ M(de) ^—-^ = [ M{de) = 1 , 

e K Ju^T J R (e-E) 2 + w J R 

the inequality (TSUI) follows upon integrating (fl9~j) . □ 

Proof of Theorem [TJ According to Proposition [1] it remains to bound EM|' u ' Q (£ ,£i) for 
a > 1. We further split this into two contributions: 

f El dE i 
M T (1) = ]T |n|W ^E|Gj +T (0,n)| 

0<|n|<(logT)£ ^• fc0 
11 



2 



and My (2) corresponding to the sum over (logT) 13 < \n\ < T a . The bound ([211) with N = 
(logT)^ immediately gives Mr(l) < (logT)^ 9 . The second contribution can be bounded using 

m 



64 r 4 T 5 

n H e 

7T 

n>(logT)/3 



m 7 (2) < ¥I_L y — 



which is bounded by a constant C(/?, g) due to Lemma [2] as long as (3 > 2. □ 

This above proof applies equally well to the half-line problem with arbitrary boundary 
condition, the only difference being that the IDS Af has to be replaced by the 'E(fi UJ (dE)) where 
/i w is the spectral measure of |0) and H. 

5 Probabilistic estimates near a critical energy 

In this section, we first derive more quantitative versions of Lemmata [3] and H] by replacing the 
input ([()]) and ([7]) by the estimates (TTTT) and (TT2"j) . However, these estimates are not sufficient 
for the proof of Theorem [3J In fact, one can further improve the Lemmata by replacing the 
uniform upper bound ffTTl) by a probabilistic one, deduced from a large deviation estimate from 
[JSSJ recalled below and showing that the transfer matrices grow no more than given by the 
Lyapunov exponent with high probability. For sake of notational simplicity, we suppose that 
E c = 0. Furthermore, according to (fT2"|) and the Thouless formula we may choose positive 
d < D and e such that 

j(z) > de 2 , for z = e + 15 with |e| < e . (21) 

In order to further simplify notation, we also assume that 5, e > even though all estimates 
hold with \S\ and |e|. 

Lemma 6 For z = e + 15 with 5 < e < eo introduce the set 

Q N (z) = {ujeQ J ||TJ(iV,0)|| 2 > e d ^ N } . (22) 

Then 

P(Q N (z)) > Cl e . 

Proof. This is an immediate corollary of Lemma [3] with c\ = d/(2co), where cq is introduced 
in <m. □ 



e J N 



Lemma 7 Let z = e + 15 with 5 < e < e . Then there exists a constant c 2 such that the set 

Q N (z) = \ uj e Q max ||7^(n,0)f > e 

I 0<n<N 

satisfies 

P{£l N {z)) > 1 - e~ C2N . 
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Proof. Let us split N into pieces of length No and follow the proof of Lemma H] invoking 
Lemma [6] instead of Lemma [3f giving 



uj e Q 



max || T* (n, 0)|| 2 > e 



■d^ N, 



o 



0<n<N 



> 1 - (1 - ae) N o > 1 



□ 



Choosing N = ^ shows that one may take Ci = c\d/2. 

Certainly other choices of N are possible in the previous proof, but the present one leading 
to Lemma [7] implies the following estimate, which is sufficient in order to deal with one of the 
terms in the next section (a boundary term of energies close to eo). 

Corollary 1 Let z = e + i5 with 5 < e < e . There is a constant C3 such that 

1 



E 



< e 



max <| n |<Ar || T z (n,0) || 2 
Now we turn to the refined statements and start by recalling the following 

Theorem 4 |JSSj Suppose that E(e 2M?CT ) 7^ 1 and < a < |. Then there exist constants 
C4,c 5 ,Cq, C7 such that the set 



"TV 



\ cueQ max \\T* +t5 (n,m)\\ < e C4 V \6\ < c^N' 1 , |e| < N~^~ a 1 , 

I 0<n,m<N I 



satisfies 



P(O^) > 1- c 6 e 



-c 7 N a 



For a fixed energy z, this can be extended to length scales N beyond the localization length 
(inverse Lyapunov exponent). 



Lemma 8 For z = e + i5 with 5 < c$e 2 and e > N 2 a j the set 



u e Q 



max \\TJ(n,m)\\ < e 

0<n,m<N 



CANf 2(l-2 a ) 



satisfies 



P(tt%(z)) > 1- c 6 Nee 



CT€ 



Proof. We split N into pieces of length N = e 2a + 1 . The condition e > N 2 a insures 
that iVo < N. By the stationarity, on each piece we may apply Theorem H] (with N instead of 

N) because e < N 2 ° '. For the jth piece, denote the set appearing in Theorem H] by Q^} 3 . For 
any ui G fl . , n_^n 3 one then has for 5 < c^Nn 1 (and hence also 5 < c 5 e 2 ) the estimate 

max \\TJ(n,m)\\ < e C4 ^ < e ^ N - 2(1 ' 2a) . 

0<n,m<N 
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Therefore D,-_i jv Qlf C Q%(z). We hence deduce from Theorem 0] that 

i — ■ L i — l M JV 

P(^(^) C ) < c 6 ^e- C77V o < c ,Nee- c ^-\ 

iV 

which is precisely the statement of the lemma. □ 

This last lemma can now be used in order to improve Lemma [6] in the range 5 < ce 2 . 

Lemma 9 Let z = e + i5 with 5 < c^e 2 and N~^~ a < e < N~ a , < a < \. Then the set 
Qn{ z ) defined in (|22|) satisfies for some constant c§ 

P{n N {z)) > c 8 e Aa . 

Proof. We argue as in the proof of Lemma [3j Let us set again P = P(fijv(z)), and estimate 
separately the contribution from the complement of Qn(z), ^n( z ) an d it's complement. Due 
to ( l2~TI) . Lemma [Eland the a priori bound (jlTj) (used on the complement of Q^(z)), it follows 
that 

de 2 < (l-P)lde 2 + Pc 4 e 2(1 - 2Q) + 2 c c 6 N e 2 e" C7 ^ . 



2 

Hence 



4Q d-Ac oC6 Ne- c ^' 
P > e - 



2c 4 - de Aa 

The hypothesis e < N~ a implies the result (it would be enough to assume e < log(iV) p for some 
p). □ 

As the final preparatory step for the next section, we improve Lemma [7] in the range 5 < ce 2 , 
by invoking Lemma [9] in its proof. 

Lemma 10 Let z = e + i5 with 5 < c 5 e 2 and N~^~ a < e < N~ a , < a < \. Then the set 



tt%(z) = { ,c e o 
satisfies for some constant eg 



max ||TJ(n,0)f > e 

0<n<N 



N l-a e 2(l + 2a) 



P{tt%{z)) > 1 - e 



-eg iV" 



Proof. Splitting N into pieces of length N and arguing exactly as in Lemma H] invoking 
Lemma instead of Lemma [31 we obtain 



uj e Q 



max^llTJ^O)!! 2 > e^ d ^ N ° \)>l- (l-c 8 e 4 *^ 



Choosing N = (2N 1 a e ia )/d shows that one may take eg = c$d/2. □ 

Lemma [10] implies the following estimate, which is the main result of this section and will 
be used in the next one. 

Corollary 2 Let z = e + %5. If 5 < c 5 e 2 and N~^~ a < e < N~ a , < a < \, one has 

max <| n |<iv || T z (n,0) \\ J 
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6 Proof of upper bound for random polymer models 



In this section we complete the proof of Theorem [3l For this purpose, we follow the strategy 
discussed in Section |3] and consider M^' a °' ai (Eq, E±) defined as in ( fl6l) . but with a disorder 
average E. By Proposition CD it is sufficient to bound M!^ 0,1+a (E , E\) for a > if (E ,Ei) 
contains the spectrum. Moreover, energies bounded away from critical energies have a strictly 
positive Lyapunov exponent |BG| . By the results of Section HI these energies hence lead at most 
to logarithmic growth in time, and therefore give no contribution to the diffusion exponents 
(3 q . Thus we are left to deal with energy intervals around the critical energies. All of them 
are treated the same way, so we focus on one of them. We suppose that E c = and consider 
only the energy interval [0, e ] with e chosen as in (|2T]) ; the other side [— e , 0] is again treated 
similarly. Furthermore we split the contribution as follows 

il4'°' 1+a (0,e ) = My 0,1+a (0, T~ v ) + Mf' 1+a {T-\T- a ) + Mf' 1+a {T- a ,e ), 

where r\ = min{g, -}. This is a good choice due to the following lemma, showing that the 
contribution M|,' 0,1+a (0, T~ v ) is bounded by the diffusion exponent as given in Theorem [3) 

Lemma 11 For some constant C\, one has 

M|'°' 1+a (0, T-v) < d T q ~ v+aq . 

Proof. By ( TT9|) we have 

M |Ai+"( 0;T -?7) < Tq (i+a) / de^mB^e + iT- 1 ) . 

Jo 

The estimate now follows from (fT5|) and Proposition [3] in the appendix, with eo = T~~ v . □ 
Next let us consider the boundary term M^°' 1+a (T~ a , eo). 

Lemma 12 For some constant C2 = C*2(a), one has 

M«' 0,1+Q (rVo) < T 4ag + C 2 . 

Proof. Let us split the sum over n appearing in the definition of Mf' 1+a (T~ a , e ) into one 
over \n\ < T Aa and the other over \n\ > T 4a . The first one can be bounded by T 4aq using (120]) . 
In the second one we apply Proposition [2] combined with Corollary [1] in order to bound the 
Green's function. This shows M^ ,0,1+a (T^ a , eo) is bounded above by 




In the second term, the sum over n is bounded by Lemma [2j As e > T~ a , it follows that the 
second term is bounded by CT p e _C3T< * for some C, p > 0. Hence this term gives the second 
contribution in the bound. □ 
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Lemma 13 For 77 = min{g, -} and some constants C3 = C^(a) and C4, one has 

r«-3 g > I , 



Mf'°°(T~ r ',T- a ) < C 3 + C 4 T 69a 



r a(2 9 -l) g < 1 



Proof. We first split M^' 0,co {T~ v , T~ a ) into two contributions M T (1) and M T (2), the first 
containing all the summands with \n\ smaller than the (energy dependent) localization length: 



M T (1) = T % £ M«E|G^(0,n)| 2 , 

1/7 l<\n\<e~ 2 T ia 

and the second My (2) containing the sum over |n| > e~ 2 T 4a corresponding to the summands 
beyond the localization length. M-r(l) is bounded using ffTUl) : 

de_ e - 2q 1 



T (E — e) 2 + T- 



M T (1) < T 4Q9 y N{dE) J 
In order to bound the factor e~ 2q , let us split the integral over e into pieces: 

77 — a 

M T (1) < T 4aq V T 2q{ ^- 1)a) / J\f(dE) / - 



77 — a 



T (£-e) 2 + T" 



Using Proposition [3] we obtain 

77 — a T) — OL 

OL OL 



2q)aj 



For g > 1/2, we use the bound T^~ 2q ^ a ^ < 1 showing that the sum is bounded by 2 = ^. 
For g < 1/2 we bound each summand by T^ 2q ^ a ^ < T^ 2q ^' n ^ a K This gives the second 
contribution in the lemma. 

It remains to show that M^(2) < C3. Due to Proposition [2] and Corollary [2J 

M T (2) < f°- ]T |n|«16 r 4 T 6 (V^ 1 ^ M 1 "* + g-o.M-N 

^ 7rT |n|> e -»T*«» 

Using Lemma [2] it is now elementary to bound Mt(2) by a constant. □ 

Combining Lemmata [TTj [12] and EH and recalling that a can be taken arbitrary close to 0, 
proves Theorem [31 
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Appendix: estimates on the Borel transform 

In Section [6] we used well-known estimates on the Borel transform B^(z) = J fi(dE) (z — E)~ x 
of a measure \i. For sake of completeness we provide a short proof. 

Proposition 3 If a measure \i satisfies at some E the bound /jl([E — e, E + e]) < Ce for all 
e > 0, then for any finite positive 5 and to 

SsmB^E + iS) < - C, / de^mB^E + e + i5) < tt 2 C e . 

Proof. One has, uniformly in S, 



^smB^{E + i5) = I fi(de 



e — E) 2 + 5 2 











eel 


\e-E\< 









t 







<: C) I ^ dt \l- -5 2 < C ! dx x/7 ' 



t Jo + 



which gives the first inequality. For the second one, let us bound the indicator function 



X[_ eo>eo ](e) above by e2 ^° a . Then one can use the stability of Cauchy distribution to obtain 

[ ° deSmB^E + e + 15) < [ de ^ C ° 2 f 
Jo J e + e o J 



2ne J jj,{de) 



(E + e - e) 2 + 6 2 



(E-e) 2 + (6 + e ) 2 ' 
= 2ire $mB ti (E + i(6 + e )) , 

and thus the second inequality follows from the first one. □ 
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